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Abstract 

We show that metastable states in Winter model can be related to the eigenstates 
of a particle in a box by means of renormalization and mixing. 
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1 Introduction 



The theory of unstables states in quantum mechanics [H [2|, [3l HI O [6l [7] has apphcations 
virtually in any branch of physics: statistical and condensed matter physics [H [9l [TOl IHl |12] , 
atomic and molecular physics [Tl] , nucler physics [12] , quantum field theory and particle 
physics [16j and so on. 

In this note we show that the particle states inside the cavity (see eq. ([5])) of the Winter 
model [ITJ [IHl [m |20l |2T] can be obtained from the states of a particle in a box (see eq. ([3])) 
by means of renormalization and mixing. Winter model describes the coupling of a cavity 
with the outside and is given, in the notation [21], by the Hamiltonian 

in the half-line < x < oo and with vanishing boundary condition in zero: iIj{x = 0,t) = 0. 
6{y) is the Dirac (5-function and g is a positive coupling. The metastable states are nothing 
but wave packets initially (i.e. at t = 0) concentrated inside the cavity (0 < a; < tt, 
see eq. (jl])). The time evolution of the resonances is controlled by wave propagation and 
imperfect multiple reflections on the cavity walls, leading to a leakage of the wave amplitude 
outside. For small coupling, < g 1, there is weak coupling of the cavity with the outside 
and resonant long-lived states come into play. The idea is that, by means of them, we can 
describe the dynamics of the particles initially inside the cavity as if the outside did not 
exist. Let us stress however the difference between the "true", exact, eigenstates and the 
resonance states of ([1]): the former have real energies and lie in the continuum spectrum, 
while the latter are normalizable states with complex energy and describe dynamics in a 
simple but approximate way for a finite amount of time only p[^. There is however a 
remarkable connection between these two set of states: the continuum eigenfunctions show 
for quantum numbers G M close to the resonance ones k^^\g) G C marked peaks inside 
the cavity compared to the outside [181 121] . 

Our results involve a conjecture based on a second-oder computation in g extending the 
0{g) results in [21], in which we repeated the original Winter's computation, finding non- 
diagonal contributions in the time evolution of the resonances which were absent in [T7j (see 
eq. ( 123|) in the next section). These terms have a small strength 0{g) <C 1 compared to the 
diagonal ones, but decay slowly in time, with the smallest decay width F*^^) (see eq. ( 12T|) ). 
These contributions therefore dominate the evolution of all the excited resonances (/ > 1) 
at large times and cannot be neglected. In this note we show that such contributions can in 
principle be "rotated away" by means of a unitary transformation which accounts for most 
of the dynamics. The latter is actually one of the renormalization effects occurring when 
we switch on the coupling of the cavity with the outside (i.e. we go from (7 = g ^ Q). 

To show the relevance of such interpretation, let us briefly discuss the importance of 
renormalization in understanding the dynamics of many physical theories. It is common 
practise in physics to compare a given physical system with a simplified one in which some 
interactions are omitted. These interactions can be related to control parameters in ex- 
periments involving, for example, external electric or magnetic fields, or can be treated 
theoretically as variable quantities. The idea of renormalization is that switching on an 
interaction in a physical system has, as one of the main dynamical effects, that of modify- 
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ing the parameters of the starting, non-interacting, system; once renormahzation has been 
made, the residual effects of the interaction are substantially weaker than before renormal- 
ization. When the coupling setting the strength of the interaction gets too large, there is 
usually little connection between the free and the interacting system and renormahzation 
often looses its meaning. 

In condensed matter physics, renormahzation is related to the so called "adiabatic con- 
tinuity principle" |T2]: by adiabatically (i.e. slowly) turning on an interaction, the free 
states of the system go 1 — 1 onto the interacting states by means of a flow of the param- 
eters such as masses, couplings, etc. A typical example is the normal Fermi liquid, i.e. a 
system with a repulsive interaction among electrons in the Fermi sphere 0. In quantum field 
theory ESI [211 125], the relation between free parameters and interacting ones (masses, 
couplings and field normalizations) is often singular because of the lack of intrinsic energy 
scales cutting off the quantum fluctuations at large energies: that leads to the well-known 
ultraviolet infinities. 

In classical physics, renormahzation is usually implemented by the method of "multiple 
scales" In the case of a free anharmonic oscillator, for example, renormahzation 

amounts to the "absorption" of secular terms into a shift of the harmonic frequency. These 
terms are formally resonances produced by forcing terms occurring in the perturbative 
expansion, are incompatible with energy conservation and spoil the convergence of the 
perturbative expansion at large times. After renormahzation, such strong-coupling effects 
completely disappear, only a small coupling between the harmonics is left and a uniform 
approximation in time is obtained. 

Let us remark that the adiabatic continuity principle is subjected to relevant violations. 
Let us quote, for example, the cases of the energy gap in the Bardeen-Cooper-Schriffer 
(BCS) theory of classical superconductors [28] or the mass gap in massless Quantum- 
Chomo-Dynamics (QCD) j2^. These phenomena are tipically characterized by functions 
which have an essential singularity when the interaction coupling g goes to zero, of the form 
gi/s < 0, making non smooth the connection between the interacting system and the 
related free one. In these cases, the relation between the free system and the interacting 
one is highly non trivial and the residual interaction is of nonperturbative character. 

Let us end the introduction by observing that, even though renormahzation is imple- 
mented and interpreted in quite different ways in different contexts, it is an ubiquitous 
phenomenon in physics — like the unstable states cited above — a thing which certainly 
could not be expected a priori. 



2 Second-Order Computation 0{g'^) 

The eigenf unctions of a particle in a box of length L = vr with Hamiltonian 



See footnote m 
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are given, as well known, by 



4')(x,t) = J^sin(/x)e-^''*, (3) 



where / = 1, 2, 3, ■ • ■ is a positive integer and < a; < tt. We study the time evolution of 
wavefunctions il)^^\x,t) which coincide at t = with the free eigenf unctions in eq. ([3]) in 
the interval x G [0, vr] (the cavity) and vanish outside it: 



,j,{i)(r {]\ - / v2/vrsin(/x) for < x < tt; , . 

^ ^^'^^ ~ \ for7r<x<oo. ^ > 

In a time slice between a pre-exponential small-t region [131 [13 HI] and a post-exponential 
one related to a power-hke decajl^ [171 120] , these wavefunctions describe mestastable states 
for < 5^ -C 1 evolving as the following sum of exponentials: 

oo 

^%{x,t-g) = Y,U{g)U^-\x,t-g). (5) 

ra=l 

Eq. (|5|) does not reproduce indeed the initial value in eq. ([1]) for t = 0. It is clear that 
non-exponential contributions do not have a resonance interpretation: they constitute an 
intrinsic limit of the scheme. 

The infinite mixing matrix U{g) has the following exact expression: 

^^^^^ ^(-ir+^2/fc(")(^) [l-27:igk^-){g)\"^ 1 

[P - m {gf]\l + g[l- 2m m {g)] } ^^"^ (^) ' 

where we have chosen the branch of the complex square root with 1^/^ = 1. The function 
k^'^\g) is the branch with A;*^"')(0) = n of the multi- valued analytic function k{g) satisfying 
the trascendental equation [2l] 



g)='- + ^ [exp {2mk) - 1] = 0, (7) 

i.e. satisfying b{k{g),g) = 0; the expansion of k^^\g) for small g is given in eq. f ITI?]) . Let us 
remark that eq. has been obtained by using eq. d?]) to eliminate the exponential terms. 
The function ^^"'\x,t;g) is the n-th pole wavefunction: 



^ ^Z^-'> {g) sin {g) x] E^^^ {t;g). (8) 



We have multiplied and divided by the renormalization constants Z'^"'\g)^s which, contrary 
to [21j, are still undefined and depend in general on n: 

Z(")(^?) = l + a(")^? + 6("V + 0(^?3), (9) 



The explicit expression for the bounds can be found in [171 HI] • 
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where a'-"'' and a^""^ are complex coefficients to be determined. The pole wave-functions 
evolve diagonally in time with the factors 



(10) 



Note that E^'^\0;g) = 1, as it should. Since the energies are complex for g ^ 0, it is 
convenient to split them into real and imaginary parts as: 



-^''Kg) = (k^'-Hg))' = u^'^Kg) - '-T^'^Kg), 



where u^"'\g) is the frequency and T^'^\g) is the decay width of the pole state n. 

To simplify formulae, it is convenient to introduce matrix notation. Ley us define an infinite 

column vector containing all the pole states 



E{x,t;g) 



( i^Kx,t-g) \ 
e^Kx,t-g) 



e-\x,t-g) 



\ 



(12) 



/ 



and an infinite diagonal matrix representing the evolution of the pole states 

8{t-g) ^ dmg{E'^'\t;g),E'~'\t;g),--- ,E^-\t;g),---). 
In more standard notation: 



(13) 



£it;g) 



f E('){t;g) 

E(^\t;g) 



\ 



E^^\t-, 



\ 





(14) 





The temporal evolution of the pole states can be rewritten in matrix notation as: 

'^{.x,t;g) = £{t;g) H(x,0;5() . 



(15) 



Similarly, let us define an infinite column vector containing the metastable wavefunctionsjl 



^{x,t-g) 



( iJ^^\x,t;g) \ 
ij^^Hx,t;g) 



iP^^\x,t-g) 



\ 



(16) 



^ In general, let us denote the vectors with uppercase greek letters and their components with the 
corresponding lowercase letters. 
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For t = one obtains a list of the initial conditions for all / = 1, 2, 3, 

/ sin(a;) \ 
sin(2x) 

^(x,0) 



[see eq. 



sm nx 



v 



(17) 



Eq. ([5]) reads in new notation: 
By inserting the small-^f expansion for the poles, 

one obtains the following explicit expressions for the coefficients entering the above wave- 
functions: 

T^-\g) = 4vrnV[l + O [g') . 

Similarly, the mixing matrix reads: 



■exp{,x,t;g) = U{g) E{x,t;g) = U{g) £{t;g) E{x,0;g). 



k^'^\g) = n-ng+ {n- inn^) g^ + (^-vr^n^ + Sinn^ -njg^ + O (g^) 



(18) 
(19) 



(20) 
(21) 



U{g)ln = Sin 



In 



l--g- 
2^ 

-1)"+' 2nl 



7r%2 .3, 1,2 
6 2 4 ' ^ 



1 



J Z^%g) 



+ 



(22) 
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{l-g + irmg) + g'^ 



'_lY+i+i2nl (/2 + 



n 



1 



z(-){gy 



where we have separated the diagonal (n = I) out of the off-diagonal (n ^ I) terms and 6in 
is the Kronecker S. As discussed in the introduction, the diagonal term 0{1) on the r.h.s. 
of eq. ( 12^ has been originally derived in [T7], while the 0{g) non-diagonal contributions, 
involving the real antisimmetric matrix 

(1 - Sin) (-l)'+"7#^, (23) 



(1) 

In 



P 



have been obtained for the first time, as fas as we know, in [21] . Let us end this section by 
noting a typo in the statement after eq. (19) of [2T]: the resonance states decay exponentially 
for t +00 but grow exponentially for x — )■ +oo — the well-known Gamow phenomenon 



Exponentiation Conjecture 

Let us now compare the results above, which are exact to 0{g'^) included, with an expo- 
nentiation conjecture. By using the identities 



l,oo ^ 
U2 _ 



l.oo 



J^2 _ jyip. 



Am? 



and ■ „ . n • 

t-'ik^-mPf 12 16m2' 



(24) 



^ These identities are obtained from Ylik^in"^ ^/(^ — m) — and J^k^m'^^ VC*^ ^ ™)^ = 7r^/3 respec- 
tively, by splitting the sums into positive and negative indices and rearranging them in order to have a 
single sum. The first identity can be found in [50] . 
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which hold for m a positive integer, one obtains 

2 i^' ')ln = (1 - ^'") (-1) (/2_^2)2 - ^ + g J " ^^S) 

Note that the squared matrix is symmetric, as it should, being the square of an antisim- 
metric matrix. We now assume that to 0{g^) 

Uig)=exp[Kig)]+g'H('\ (26) 

where 

K{g) = gK^'^ + g^K^'^ (27) 

is a real antisimmetric (therefore antihermitian) matrix and H^'^^ is a "remainder" ma- 
trix, taking into account residual effects not included by the exponentiation of K{g). The 
equations (!26l) turn out to be compatible and the results read: 

Z^-\9) = l-'-+(l-^nln]g^ + 0{g')■, (28) 



2r)/ 

K{g)i^ = ^(l-5,„)(-l)'+"^^A_ + . (29) 



P — n 



2 

2/n2 



f + 0{g'] 

g'Hi^ = i^g\l-5,^){-ir--^^. (30) 



P —n 

Note that if^^) is a purely imaginary matrix, which is neither simmetric nor antisimmetric. 
As reasonable in the description of quasi-particle states, decay effects are confined only 
approximately in the widths V^'^\g). 

We now generalize eq. (|26|) to all orders in g conjecturing that 

U{g) = eMK{9)] + H{g), (31) 

where K{g) is an antihermitian matrix possessing a series expansion in g starting in first 
order, 

oo 

K{g) = Y,9'K^'\ (32) 

s=l 

while H{g) is an infinite matrix starting in second order in g and collecting all the structure- 
less terms, 

oo 

H{g) = Y,9'H^'^. (33) 

s=2 

We have explicitly computed K^^\ f^^^^ and H^'^\ We conjecture that eq. flHT]) provides a 
convenient organization of the series in g of U{g), i.e. that it resums a large part of its 
series. To be more exphcit, given as an input a truncated expansion of U{g), 

N 

U{g) = Y,9''U^''\ (34) 

n=l 
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where < oo is the order, we beheve that the expression flHTl) . involving a reorganization 
of the terms in eq. (IMI) . is a much better approximation to the exact result than eq. (IMI) 
itself. 

It is interesting to compare the Z's in eq. (128|) with those obtained by imposing that the 
pole states fll2l) are normalized to one. By squaring the pole wave-function and integrating 
over X G [0, vr], one obtains: 

Zi:^M = l-l+(j-^ln^9' + 0{9'). (35) 

The above renormalization constants coincide with the previous ones (1281) to first order in 
g, while the imaginary contributions appearing in second order are different. 
It is no surprise that the terms g'^H^'^^ + g^H^^^ + ■ ■ ■ without an exponential structure 
come in second order, where they are purely imaginary: they represent effects related to 
the decay width, which is 0{g'^), and the latter is not a renormalization effect (see later 
for discussion). The terms in eq. fl55]) are the analogue of the imaginary contribution to 
ZW(^) on the r.h.s. of eq. (|28D. 

Finally, let us note that the appearence of the combination 

9 (l - f ) (36) 

in the first line in eq. ( l22|) . as well as in the expression for K{g), seems to suggest a (finite) 
renormalization of the coupling g of Winter model. 



3 Discussion 

In order to express the metastable wavefunctions of the Winter model in terms of the eigen- 
functions of the particle in a box, one has to diagonalize the time evolution. That is achieved 
by "counter-rotating" the vector containing the initial conditions, i.e. by considering the 
evolution not of \E'(x,0) but of 

^x,0-g) = U-\g)^{x,0). (37) 

By using the first equality in eq. ( |T8l) . it is immediate to show that 

<l>(x, t; g) = E{x, t; g) = S{t; g) S(a;, 0; g) . (38) 

By looking at the vector equation (1371) component by component, the new initial conditions 
read 

0«(x, 0-g) = J-Y^ {U-\9)),, Mnx), (39) 

' ^ n=l 

each evolving as a single pole wavefunction: 

(P^'\x,t;g) = ^^'\x,t;g). (40) 

The "experimental meaning" of eq. (137|) or eq. (139!) is clear: in order to observe a diagonal 
time evolution as in the free case ([3]), one has to prepare the initial state as the coherent 
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superposition of free eigenfunctions given by eq. fl57|) or by eq. f l5^ . In the case of excited 
states, / > 1, the superposition in eq. ( 137|) or (139!) has also the effect of subtracting the 
contributions from smaller /'s, which decay slower in time and therefore tend to dominate 
the evolution, as discussed in [21]. If the matrix U{g) is computed in an approximate way 
(as a truncated power series in g for example), there is a small residual contamination 
in the time evolution of the l-th state from the lower ones, which becomes substantial 
asymptotically in time. In other words, the problem of isolating the l-th mode for all times 
can in principle be solved only with an exact computation of U{g). 

Let us now discuss the renormalized wavefunctions (f)^^\x^t] g). The main qualitative 
difference between the free case and the interacting one is that in the latter case there 
are non-zero widths. The appearence of an imaginary part in the ah initio real energy is 
a second order effect in Once a non-zero width is allowed, the key point is that the 
(j)^''\x, t; gYs have a similar form to the eigenfunctions of the free system iPq{x, t) in eq. ([3]). 
The differences between the free case and the interacting one, so long as < <C 1, can be 
relegated to small modifications of the parameters entering the free wavefunctions iI)q\x, t). 
In other words, switching on the interaction, i.e. going from 5f = OtoO<(7<Cl, produces 
finite renormalizations only. Let us discuss these renormalizations in turn: 

1. the normalization coefficient Z^^\g) has a modulus smaller than one for < (7 <C 1 
and reduces to 1 in the free case (7 = 0; it has a first-order correction in g and is 
the analog of the field renormalization constant Z in quantum field theory, which is 
in general smaller than one [2T]. Unlike the most common cases (consider QED for 
example), Z^^\g) is not real because the one-particle states are unstable; 

2. the wave- vector k'^"'\g) is renormalized to first order in g by the interaction and 
reduces to the free case for ^ ^ 0: A;«(0) = /. It acquires an imaginary part in 
second order in g, related to the decay width. That implies the disappearence of 
the node of the wavefunction around x = n and a (small) exponential growth of 
(f)^^\x^t]g) by going from the impermeable wall in x = toward the permeable one 
in X = tt; 

3. the real part of the energy u^''\g) is also renormalized to first order in g by the 
interaction and reduces to the free case for g ^ 0: a;*^'^(0) = P. Note that the free 
dispersion relation u = k"^ is not renormalized in first order. 

Let us make a few remarks. 

• We do not expand in powers of g the wavefunctions 0*-'^(x, t; g), but only the parame- 
ters k^'^\g), u^'^\g), etc. entering them through the functions appearing in il)Q\x,t). 
That implies that we are resumming classes of higher order corrections in g in the 
wave function, in the spirit of renormalization in quantum field theory [221 1231 [211 125] 
and statistical mechanics [3T] or the method of multiple scales in classical physics [2Z] ; 

• The decay widths grow faster with increasing n than the frequencies: 

^(")(^)ocr^^ u^''\g)o^n^. (41) 

^ Non zero widths are clearly not in contradiction with the unitarity of the fundamental thery because 
we are looking at a subsystem, an "open" system [H [TH] . 
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Since our renormalized theory has meaning only for 

T^^\g)<^J-\g), (42) 

we cannot take n too large. Therefore, while in principle the state vectors and the 
evolution/mixing matrices are infinite, in pratice for any fixed g one has to make a 
truncation in n according to the condition (1421) [§. This limitation is also reasonable 
from physics viewpoint: high energy particles pass through the barrier in a; = vr 
without difficulty and therefore there is no sense in including them to describe the 
dynamics inside the cavity. By restricting on n one is also cutting off small wavelenghts 
\<2t{ /n and therefore is limiting space resolution. 

4 Conclusions 

We have shown that the resonant states of the Winter model can be related to the states 
of a particle inside an impenetrable box by means of renormalization of the parameters 
entering the "free" eingenf unctions and mixing of the states, after allowing for non-zero 
widths. The less trivial aspect of this renormalization procedure is related to a term which 
was overlooked in [17] and has been found in [21]: the initial resonant states (at t = 0) 
have to be properly "counter-rotated" by means of the infinite matrix U{g)~^, the inverse 
of the state mixing matrix U{g), in order to have a diagonal, exponential time evolution. 
We have explicitly computed U{g) to second order in ^ 1 and we have conjectured a 
general form for U{g) involving a unitary matrix plus a matrix collecting residual decay 
effects (see eq. (13T1) ). The unitary matrix is the exponential of an antihermitian matrix, to 
be computed order by order in g — a typical phenomenon in quantum field theory. 
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